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ABSTRACT 

Aims. We derive simple estimates of the maximum efficiency with which matter can be ejected by the propeller mech- 
anism in disk-fed, rotating magnetic neutron stars. Some binary evolution scenarios envisage that this mechanism is 
responsible for expelling to infinity the mass inflowing at a low rate from the companion star, therefore limiting the 
total amount of mass that can be accreted by the neutron star. 

Methods. We demonstrate that, for typical neutron star parameters, a maximum of ^ pr o < 5.7(P-3) 1//3 times more 
matter than accreted can be expelled through the propeller mechanism at the expenses of the neutron star rotational 
energy (P-3 is the NS spin period in unit of 10 -3 s). Approaching this value, however, would require a great deal of 
fine tuning in the system parameters and the properties of the interaction of matter and magnetic field at the magne- 
tospheric boundary. 

Results. We conclude that some other mechanism must be invoked in order to prevent that too much mass accretes 
onto the neutron stars of some low mass X-ray binaries. 

Key words, accretion: disk - stars: close binaries -stars: neutron -pulsars: general 
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One of the possible endpoints of the evolution of the com- 
pact object in a low-mass X-ray binary (LMXB) is a mil- 
lisecond pulsar (MSP), a rapidly spinning neutron star 
(NS) with weak surface magnetic field (~ 10 8 — 10 9 G). 
Such NSs are thought to be spun-up to millisecond pe- 
riods by the accretion of mass and angular momentum 
from matter in a disk formed through Roche-lobe over- 
flow of a low-mass companion star (fo r a review see e.g., 
iBhattacharva fc van den Heuvell . ll99ll ). After the accretion 
phase has ended, the NS may turn on as a radio MSP. In 
recent years, the prediction that LMXB host NSs with spin 
periods in the ms range has been confirmed through the 
discovery of 300 - 600 Hz nearly-coherent oscillations dur- 
ing type I X-ray bursts from some 10 LMXBs of the Atoll 
class (e.g., IStrohmaver. Zhang fc Swankl Il997l ). This was 
followed by the detection of coherent X-ray pulsations in 
the persistent emission of seven low luminosity transients 
of the same class. The spin frequencies lie between 180 
and 600 Hz and the orbital periods between 40 min and 
5 hr (see review by IWiinanda . 120051 ). For the first time, 
the predicted decrease of the NS spin period during accre- 
tion was measured in the accreting MSP IG R J00291+5934 
(|Falanga et all . 120051 : iBurderi et all . l2006fh This provided 
a strong confirmation of the theory of 'recycled' pulsars in 
which the old NS in LMXBs become millisecond radio pul- 



sars through spin-up by transfer of angular momentum by 
the accreting material. 

All accreting MSPs are X-ray transient systems. They 
spend most of the time in a quiescent phase, with X-ray 
luminosities of the order of 10 31 — 10 32 erg s~ x , and oc- 
casionally show weeks to months long outbursts, reach 
X-ray luminosities of 10 36 
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which coherent X-ray pulsations are observed. The mag- 
netic fields inferred for these NS are of order 10 8 — 10 9 G 
(e.g.. ICampana et al l. Il998t iPsaltis fc Chakrabartvl 119991 : 
iDi Salvo fc Burderi I2003T T These findings added a great 
deal of confidence in the likely connection between accret- 
ing NS and millisecond radio pulsars. 

Models for the evolution of LMXBs as progeni- 
tors of radio MSPs still involve significant assumptions 
about the transfer of mass from one star to the other 
and the role of dynamically influent NS magnetic fields 
(|Podsiadlowski. Rappaport fe Pfahl l2002f ). The MSP for- 
mation depends on the amount of matter accreted from the 
companion which, in turn, is likely limited by mass ejec - 
tion from the binary. ICook. Shapiro fc Teukolskl (|1994ah . 
following a full GR treatment, show that a NS of 1.4 
M Q has to accrete ~ 0.2 — 0.5 M ^ to reach per i ods i n 
the range 0.6-1.5 ms (see also e.g.. IBurderi et all [l999). 
Most donor stars in systems hosting recycled radio pul- 
sars must have lost most of their mass during their evo- 
lution. They now appear as white or brown dwarfs o f 
mass ~ 0.15 - 0.30 M (|Taam. King fc Ritterl . l2000h . 
the prog enitors of which were likely s t ars of ~ 1.0 — 2.0 
M (e.g.. IBurderi. King fc Wyniill99d[Tauris fc SavonijeL 
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1999). Hence, a crucial issue for evolutionary models for 
MSPs is which fraction of mass lost by the companion 
(Mio S t ~ 0.7 - 1.8 M ) effectively accretes onto the NS. 
Mass ejection must be efficient, at least in some cases, to 
avoid that accretion of a large amount of mass onto the NS 
could induce collapse to a black hole. For most equations 
of state of ultra-dense matte r, this is expected to occur 
for NS masses above ~ 2 M Q (jCook, Shapiro fc TeukolskH 
1994b). Different models have been proposed to explain how 
mass can be ejected from binary system, thus limiting the 
amount o f matter that is accre t ed by the NS to about 0.2- 
0.4 (llllarinov fc Sunvaeyl. Il975t IWang fc Robertson! . 
119851 : iRomanova et all 120051 : lUstvugova et all l2006f ). 

In this work we use simple arguments to derive an upper 
limit on the efficiency of the propeller mass-ejection mech- 
anism. In section [2] we consider accretion and ejection of 
infinitesimally small amounts of matter, the case relevant 
to NS X-ray transient. Since in this case the NS parameters 
do not change significantly between each accretion/ejection 
cycle an analitic expression for the propeller efficiency can 
be derived. In section [3] we extend this limit to the case of 
accretion and ejection of finite amounts of matter. 



2. Ejector-Propeller efficiency 

Matter transferred from the companion star via Roche- 
Lobe overflow possesses enough specific angular momen- 
tum that an accretion disk forms around the accreting 
NS. The disk is truncated at the magnetospheric radius, 
R m . Different calculations, carried out by adopting differ- 
ent models and/or assumptions, produce only slightly dif- 
ferent estimates of the value of R m (a factor of 2- 3 at the 
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most, e.g., L . , 
IWa^fl987f) 

The value of i? m as a function of the system parameters 
is only of marginal relevance to the present paper. Therefore 
we consider he re the value inferred fro m simple theory (see 
e.g., Eq. (4) of lLamb fc Pethicld . Il97l : that is 



R m = 3.2 x lO 8 ^ 1/7 i4o 7 M 17 2/r cm 



(1) 



where Mi is the NS mass in units of 1 M , ^ 30 is the NS 
magnetic dipole moment in units of 10 30 G cm 3 , and Myi 
is the mass accretion rate in units of 10 17 g s _1 . 

Disk accretion onto a spinning, magnetic NS is expected 
to be centrifugally inhibited, once the magnetospheric ra- 
dius is larger than the cor otation radius, -R C oi and insi de the 
light cylinder radius, R\ c (jlllarinov fc Sunvaevl . Il975l ). The 
corotation radius, R col is the radius at which the Keplerian 
frequency of the orbiting matter is equal to the NS spin 
frequency: R co = (GAt/Ofi) 1 / 3 = 1.5 x lO 8 i^ 3 (M/M )V3 
cm, where fio is the angular velocity and P-3 is the spin 
period of the NS in milliseconds. The light cylinder radius, 
i?i c , is the radius at which the corotating magnetic field 
lines reach the speed of light c, R\ c = c/Oq. Outside i?i c 
the magnetic field lines cannot corotate any longer and a 
radiative electromagnetic field must be generated. If an ele- 
ment of matter dM ln is accreted onto a NS from a Keplerian 
disk, part of its energy is stored as rotational energy of the 
NS; this can later be used to eject a mass dM out . We de- 
fine the propeller efficiency as the ratio of the ejected to 
accreted matter ry pro = dM oxlt /dM- m . Moreover we assume 
that only the rotational energy gained by the star during 



the previous accretion phase can be used to eject matter in 
the propeller phase. 

In order to determine the rotational energy gained by 
the NS we first consider the accretion phase. During accre- 
tion the rate of rotational energy transferred to the NS can 
be written in terms of the transfer of angular momentum, 
L, at the inner boundary of the disk: 



E ro t — QqL — 2^0-^1 



(2) 



where L = Itt + n I = Ml. Here I = /3 _1 Mi?§, s is 
the moment of inertia of the NS, and I is the specific 
angular momentum of the accreted mass. For NS mod- 
els with mass above an initial value of M ~ 1.4 M Q 
()Weaver. Zimmerman fc Wooslevl , ll978t ) and well below the 
maximum mass, values of (5 are expected in the 2.3-3.1 
range for a variety of equations of state. In our calcula- 
tions we use p = 2.5, the value for a uniform Newtonian 
star. For a disk truncated at the magnetospheric radius, 
i? m , we approximate the specific angular momentum trans- 
ferred to the NS with the Keplerian value at R m , i.e. 
I ~ l m = (GMRm) 1 ' 2 . According to detailed models of 
the disk-magnetosphere interaction this approximation is 
accurate in the limit in which R m <C -Rco, whereas it 
provides an upper limit in the case in which i? m ~ R co 
and non-material (spin-down) torques due to the interac- 
tion of the NS magne tic field with disk regions outside 
R C o, become important fGhosh fc Lambl . ri979tlWand . ll987l : 
iRappaport. Fregeau. Spruiti . 120041 ). Since our aim here is 
to derive an upper limit on the propeller efficiency, we can 
safely adopt the Keplerian value of the specific angular mo- 
mentum. Equation ^ thus writes 



dE Iot = QoiGMRn^dMin - -^.R&gfigtiMin, 



(3) 



which describes the amount of energy accreted from a disk 
ending at the magnetospheric radius. For a given angular 
velocity of the NS, the rotational energy, dE TOt , is high- 
est when the magnetospheric radius, i? m , takes the largest 
value, compatible with accretion onto the NS surface: that 
is at R m = R co . Setting R m — R co Eq. ([3]) gives 



dE IO 



GM 
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dM in - -p-iRlMZdMn 



(4) 



This is the maximum increase in E mt that an element of 
disk-accreting mass dM in can cause. 

We now consider the propeller phase. According to the 
virial theorem, the energy required to eject a disk mass 
element, dM out , from a radius R to infinity, is dE C j C = 
GMdM out /2R. In order to estimate the maximum possi- 
ble propeller efficiency we suppose that ejection takes place 
from the largest possible radius at which the propeller can 
work, such that the energy required to eject is minimized. 
That is R = R ic . 

Setting dE rot = dE D - ]c and assuming that ejection takes 
place at R\ c , we have for the maximum efficiency of the 
propeller: 



M \-i 



where we used i?Ns — 10 km. 

An alternative possibility is that there is no magne- 
tosphere in the accretion phase, and the disk extends all 
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Fig. 1. The highest possible propeller efficiency, 77 plo , the 
ratio of the ejected mass to the accreted mass, versus spin 
period. The dashed and dot-dashed lines correspond to Eq. 
© (panel a) and Eq. ((7]) (panel b), for a NS of 1.4 M Q and 
2.0 M©, respectively. 



the way to the NS surface (such that we set formally 
R m ~ i?NS, where i?NS is the NS radius). This might be 
the case in LMXBs which do not display coherent pulsa- 
tions in their persistent X-ray flux. In this case 



Q (GMR NS ) X / 2 



1 r?2 r>2 
rt NS s: 



GM/2R lc 



c(3- 



GM 



GM/R lc 
-R^ s Qo , 



(6) 



1.48 x 10 5 (M/M Q ) cm the gravita- 



with i? g = GM/c 2 
tional radius. 

For a typical NS, i?NS — 6i?o-, one has 



(7) 



a result that is similar to that in Eq. ([5]). 

The above derived limits on ?7p ro = dM out / dM- ln involve 
a great deal of fine tuning with respect to the ejection of 
matter from close to the light cylinder. Indeed, despite the 
fact that the NS magnetic field lines approach c at R\ c , it 
assumed that ejection takes place at much lower speeds, 
precisely the lowest required for matter to reach infinity. 

Eqs. (J5]) and ([7]) are plotted in Fig. [I] as a function of 
the spin period and for different values of the NS mass. In 
the following we discuss the impact of the limits in Eqs. © 
and ([7]) on evolutionary models of accreting NS's binaries 
and MSP's. 
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Fig. 2. The highest possible propeller efficiency, 77 pr o, versus 
spin period for different amounts of mass accreted onto a 
1.3 Mq NS are shown by dot-dashed lines. The solid line is 
the propeller efficiency upper limit (Eq. [5]) . 



3. Discussion 

The limits that we have derived on the mass ejection effi- 
ciency of the propeller mechanism are based on the assump- 
tion that, when accretion takes place, the highest angular 
momentum is transferred to the NS, whereas ejection oc- 
curs at the expense of the NS rotational energy from the 
largest radius compatible with the propeller mechanism. 
This minimizes the loss of rotational energy. 

Somewhat different limits are obtained depending on 
whether during the accretion phase, the disk ends at the 
magnetospheric boundary or extends all the way to the NS 
surface. The latter case might be relevant to high luminos- 
ity NS systems which do not display periodic pulsations. 
The corresponding limit on r] pro (see Eq. [JJ applies only 
over a limited range of fast NS spin periods. Indeed, if the 
magnetosphere is "squeezed" inside the star surface dur- 
ing the accretion phase (formally R m ~ i?Ns)> then, as the 
mass inflow rate decreases by up to ~4-5 orders of mag- 
nitude (as in ferred from the lum i nosity NS X-ray transient 
systems, e.g.. lCampana fc Stellal . l2000f) . the magnetosphere 
can expand up to ~20 i?NS at the maximum. This is due 
to the stiffness of the NS magnetosphere, reflected in the 
weak dependence of R m on M in Eq. (fTJ). For the NS to 
be in propeller phase the corotation radius, R co , must be 
smaller than the magnetospheric radius R m (~ 20Rns)i i n 
turn implying that the NS spin cannot be longer than ~ 30 
ms (for NS with spin period longer than this, the propeller 
efficiency limit in Eq. ([JJ would thus be inapplicable). In 
this case, from Fig. [1] (panel b), we note that the maximum 
propeller efficiency cannot be larger than ~ 4 for NS with a 
spin period of ~ 1 ms. The limit on Tfp ro derived in the case 
in which there exists a magnetosphere also during the ac- 
cretion phase (see Eq. [5J applies instead over a large range 
of spin periods and its value is ~ 5, for a spin period of ~1 
ms and increases for longer periods. 

Strictly speaking, our reasoning applies only to the ac- 
cretion/ejection of (infinitesimally) small amounts of mat- 
ter. The limits we derived are thus immediately relevant 
to NSs that often cycle between the accretion and ejec- 
tion regimes. These include transient X-ray systems which 
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alternate active periods, lasting up to years, to quiescent 
intervals of up to ~ 100 yr. During each outburst cycle the 
NS paramaters (e.g., mass, spin period, magnetic dipole 
moment) change only by very small amounts. These sys- 
tems likely evolve through a very long sequence of accre- 
tion/ejection cycles, while within each cycle the efficiency 
of the propeller ejection remains limited by the values we 
have derived. It is easy to see that in this case the overall 
efficiency must be limited by the maximum possible value 
°f ?7pro) given by Eq. ([5]) independent of whether accretion 
occurs through the magnetosphere or not. 

Therefore, flunging away matter with higher efficiency 
than permitted by the above limit would require a different 
ejecti on mechanism such as the r adio pulsar ejection mecha- 
nism (|Illarinov fe Sunvaevl Il975t ) . This mechanism can op- 
erate at much larger radii (comparable to the NS Roche- 
lobe radius) where the gravitational potential is much shal- 
lower and thus the energy required to expel ma tter to in- 
finity is much smaller (e.g.. iBurderi et aTTl2001l ). 

In the case of accretion and ejection of finite amounts 
of matter, an analitical expression of ?7 pro cannot be found. 
In fact, in this case, one has to consider the variation of 
the stellar parameters during the evolution of the system. 
The equations relating the rotational energy acquired dur- 
ing accretion to the ejection energy lost during the propeller 
phase depend on the instantaneous value of the stellar mass 
and angular velocity, and their evolution can be followed 
numerically. To obtain an estimate of the maximum pro- 
peller efficiency we again assumed that the accretion phase 
takes place from the corotation radius, the ejection phase 
from R\ c and that all the rotational energy acquired during 
accretion is used to eject matter from the system. In our 
numerical integrations we neglected the decrease in the NS 
radius which results from increasing its mass through accre- 

— 1/3 

tion. In any case, using the simple relation Rns oc M ns 
we checked that this is a small effect and the induced vari- 
ations on r] plo are less that ~ 2%. The parameters that we 
used in our calculations are thos e in the evolutionary mod - 
els of lErgma fc Sarnal (|2003t) and lTauris fc Savoniiel (|1999l) . 
In the first paper the authors assume that propeller ejec- 
tion takes place from an increasing radius, as the mass in- 
flow rate toward the N S decreases. Using their parameters 
(|Ergma fc Sarnal . l2003f) . i.e. NS mass of ~ 1.4 M Q , initial 
spin period ~ 0.03 s, ~ 0.25 M Q accreted onto the NS, 
and a final spin period of ~ 3.6 ms, our calculations show 
that the efficiency is bound by the maximum value of < 5. 
This is a resonable result since accretion causes fio to in- 
crease and consequently, ejection takes place (mainly) at 
smaller values of R\ c . In this case the energy required to 
eject matter is higher and the propeller less efficient. An 
additional effect comes from the fact that the corotation ra- 
dius slightly increases during accretion due to the building 
up of the NS mass. However, our calculations show that this 
effect is small (due to the limited range in NS mass). Infact 
in all cases we obtained a lower efficiency than the maxi- 
mum value of ry pro of Eq. ([5]) , which holds for short accre- 
tion/ejection cycles. Fig.[5]shows that the upper limit of Eq. 
(O is approached only when the accreted mass on the NS is 
very small. We note that the prope ller efficiency required in 
the work of lErgma fc Sarnal (|2003f ) is ~ 2.6, consistent with 
the limit we have derived. It is more difficult to carry ou t 
a comparison with the work of iTauris fc Savoniiel (| 19991) . 
We refer to their models for relatively close binary systems 



with light donors (P orb < 10 d, M donor < 1.3 M s ). This 
is because, in those cases, the mass-transfer rate remains 
sub-Eddington during the entire binary system evolution, 
preventing the formation of super-Eddington strong out- 
flows and causing the binary system to eject matter mainly 
via the propeller effect. We use a NS mass of ~1.3 M Q , 
~0. 1-0.4 M Q accreted onto the NS and an initial NS spin 
period in the range ~ 0.001-0.006 s. Again we find the re- 
quired propeller efficiency in Tauris & Savonije (in some 
cases 77 pro < 1) is smaller than our numerically estimated 
upper limit (r? pro < 7). 

A refined treatment of the propeller eff ect for a magnetic 
accret ing ms pulsar was carried out by lUstvugova et al.l 
(2006). Using a magnetohydrodynamic approach, and con- 
sidering a detailed model for the interaction between the 
NS magnetosphere and the accretion disk, they showed 
that propeller outflows can be expected from a fast rotating 
(-Pspin ~ 1.2-1.6 ms) weakly magnetic NS (B ~ 5 x 10 s - 
2 x 10 9 G). A comparison between their results and ours 
77 pro is straightforward if one consider that all the NS ro- 
tational energy used in the propeller phase derives from 
accretion of matter at the corotation radius. Dividing their 
total ejected matter (Mj + M w ) by the above estimated 
accreted matter, we obtain ?7 pro ~ 0.4, for a 1.26 ms pulsar 
with B ~ 10 9 G. 

In summa ry the evol utionary calculat ion of 

lErgma fc SaTnal (|2003h and ITauris fe Savoniiel (|l999h 
correspond to a propeller efficiency of ~1— 3, which as such 
would be adequate to explain how the ejection of a large 
fraction of the mass transferred by the companion can 
take place in these systems. The upper limit that we have 
derived in this paper for the propeller efficiency is only a 
factor of ^2 higher. However achieving such a limit would 
require a high degree of fine tuning. This indicates that 
a propeller efficiency of order ~2-3 might be difficult to 
achieve in real systems. 

As suggested by th e deteiled propeller simulations by 
lUstvugova et al.l (|2006f ) a propeller efficiency < 1 might be 
more realistic. If this was the case, then a different mecha- 
nism to eject from the binary most of the mass transferred 
by the companion is required. The onset of the radio pulsar 
ejection mechanism is thus probably to be favored. 
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